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Abstract 

The excitation of a non-linear ion-wake by a train of ultra-relativistic electron bubble wake 
[jll-iSI is modeled. The ion-wake is shown to be a driven non-linear ion-acoustic wave in the 
form of a cylindrical ion-soliton llTT] - iT4l . The phases of the oscillating radial electric fields of the 
slowly-propagating IIZTl electron bubble is asymmetric in time and excites time-averaged inertial 
ion motion radially The electron compression in the back of the bubble sucks-m the ions and 
the space-charge within the bubble cavity expels them, driving a cylindrical ion-soliton structure 
with on-axis and bubble-edge density-spikes |5l[|6l. Once formed, the charmel-edge density-spike 
is driven radially outwards by the thermal pressure of the wake energy Ill2l . Its channel-like 
structure due to the flat-residue left behind by the propagating ion-soliton, is independent of 
the energy-source driving the bubble lUljll electron wake. We explore the use of the partially- 
filled channel formed by the cylindrical ion-soliton for a novel regime of positron acceleration 
[fl8] . OSIRIS PIC l( 25l simulations are used to study the ion-wake soliton structure, its driven 
propagation and its use for positron acceleration. 
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1. INTRODUCTION 


Plasma ions are generally assumed to be stationary in the theory of ultra-relativistic 
non-linear plasma electron waves excited as wakefields of a laser or particle beam IfH-Hbl. 
Important exceptions to this occur when the drive beam (or a witness beam being accel¬ 
erated) in a plasma wakefield accelerator becomes denser than the background plasma 
density by a factor larger than the ion to electron mass ratio ll9ll[(T0ll . Such is the case for 
the parameters of ultra-low emittance future collider designs at the TeV energy scale. Ion 
motion also becomes important in the wake train behind the initial wake oscillation or 
bubble. Understanding the long-term ion behavior |jZll|8ll is important to determine the 
state of the plasma for succeeding pulses in a future plasma collider of high repetition 
rate IfT^lfTT!! . We also show that the resulting ion-wake chaimel formed is of interest for 
positron wakefield acceleration. 

In this paper, we model the excitation of a nonlinear ion-wake in the trail of an electron 
plasma wake in the bubble regime, driven by either relativistic ally intense laser or particle 
beam energy sources. The ion-wake is shown to be a non-linear ion-acoustic wave in the 
form of a cylindrical ion-acoustic soliton. Its characteristics are similar to the solutions 
of the cylindrical Korteweg - de Vries equation (cKdV) lUTI - IlT^ . However, the ion-wake 
shown here is different from a cylindrical-KdV soliton due to several reasons: The bubble 
wake electron oscillations do not thermalize into an isothermal plasma; the ion-wake 
soliton breaks up into N-solitons as it evolves and there is an ion-density spike on the 
axis. The soliton propagates radially outwards leaving behind a flat residue resulting in 
a partially filled ion-wake chaimel. The dynamics of the ion wake is shown to have two 
distinct phases. The excitation phase where the time asymmetry of the radial electric fields 
of the bubble excite the ions into a soliton-like structure. At later times, the propagation 
phase where the soliton is driven outwards by the pressure of thermalized electrons after 
the electron oscillations undergo phase mixing ||20l. 

Representative PIC simulation results in Figj^ and Fig|^ illustrate the salient features 
of the non-linear ion-wake. Figur^ shows the excitation phase at an early time when 
the bubble wake train is executing ordered oscillations and its fields have begun to excite 
the ion-wake structure seen in Fig|T|b),(c). The radial oscillations later phase mix and 
the electrons thermalize, their thermal pressure then drives the ion-soliton outwards as 
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shown in Fig|^ As we show below, it is the longitudinal or time asymmetry of the 
radial electron wakefields that excites the ion soliton which propagates leaving behind a 
partially-filled chaimel shown in the PIC simulations in Figj^and Figj^ 

In section]^ we present considerations and assumptions made to derive the ion-wake 
model. In section |n^ using the fluid equations for ion dynamics we model the non-linear 
ion-acoustic waves as a driven cylindrical ion soliton. The driven wave models are shown 
to have two distinct phases: the excitation phase and the propagation phase. We present a 


and verify it using simulation results. In section|V]we present simulation results showing 
the dynamics of propagation of the ion-soliton driven by a radial temperature gradient. 
Finally, in section|V^we analyze the properties of a positron beam driven wakefield in an 
ion-wake chaimel. 


model for the excitation of the ion-soliton by the electron bubble radial fields in section IV 
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FIG. 1; Laser driven non-linear ion-wake at early time (t = = 0.17 f~.^) in in, = nip = 1836 plasma, (a) Electron bubble 

Wakefields in cartesian coordinates (fixed-box) with ^ = 10 driven by a matched laser pulse (vector potential tio = 4 and frequency 
oiq) with Rg - 4^. (b) Non-linear ion-wake in the form of a cylindrical ion-soliton of radius - 4^ excited behind the bubble 
electron wake in a proton plasma, (c) Transverse ion-density profile at z = 15 c/aip^. Notice that the ion density perturbation in this 
excitation phase is still building up and is a fraction of the background ion density, ^ < 1 ■ 
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FIG. 2; Electron beam-driven non-linear ion-wake at late time (t = 460 = 1.7f~}) in m, = nip = 1836 plasma, (a) 

Beam-driven ion-wake electron density in cylindrical coordinates (fixed-box). The beam parameters are nj, = 5no, Or = 0.5cjwpc, 
ffz = 1.5clcope, yi, = 38.000, these beam-plasma parameters are quite similar to (H. (b) Corresponding ion density in cylindrical 
coordinates (fixed-box). Note the N-soliton formation in the ion-density. This is seen as the ion-wake evolves further, for instance at 
2 = 60 cjcDp}. (c) Radial electron and ion density profile at z = 150 clcopg. 
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II. CONSIDERATIONS FOR MODELING THE NON-LINEAR ION-WAKE 


We recognize that to study the time evolution of a wake-excited plasma in order to 
establish the duration over which it relaxes to thermal equilibrium both collisional and 
collision-less dynamics have to be considered along with the physics of recombination 
modes such as electron-ion recombination. However, in this work the dynamics of plasma 
is modeled under the collision-less approximation. Thus, diffusion is not important dur¬ 
ing the timescales over which the ion-wake is studied. We do not discuss recombination 
except mentioning that the "afterglow" is dominated by volume recombination while 
localized effects cannot be ruled out. 

In order to formally establish the difference between the density wave processes that 
occur over collision-less timescale in contrast to the ones that start dominating under 
collisions, we show the assumptions made to arrive at the dynamics of diffusion. The 
process of diffusion is modeled with a parabolic partial differential equation which is 
deduced from the ion-fluid equations under the assumption that the inertial response of 
the ions is much faster than the collisional timescales. 

The effect of collisions is introduced as a drag force, m/n Vcoii (ti)- The collisional 
drag force modifies the ion-fluid equation of motion as, = mn + {ti)V • (vi)^ = 

±e nE - We - mini Vcou (Vi) where Vcou is the average electron-ion collision frequency and 
is obtained from the mean free path. Diffusion of plasma is thus driven by the charge- 
separation field, E and the thermal pressure, Pe while being impeded by the collisional 
drag. Upon ignoring the inertia of the ions, ^ = 0, the equation for the ion velocity 
by diffusion in an isothermal plasma is, (Vj) = ±—— E —phe characteristic 
parameter of diffusion is fhe diffusion coefficient or diffusivity D = and mobility 
p = which depend upon the collision frequency. Using the gradient of the velocity 
in the continuity equation and ignoring the mobility p, leads to a Tick's law diffusion 
equation, characteristic of a parabolic equation. 

When the mobility is retained, the fluid equation is a moment of the Fokker-Planck 
equation which is the kinetic model of the collision-driven drift and diffusion. The 
diffusion equation thus cannot support wave-like solution because such solutions are 
characteristic of a hyperbolic partial differential equation. 

The solutions of linear and non-linear diffusion equations show the evolution of den- 
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sity profile by diffusion and can be obfained using the self-similar formulation. The 
self-similar solutions show the spatial and temporal evolution of the density to be ex¬ 
ponentially decaying. In the non-linear case, the density can have a sharp-front as it 
decays. However, a soliton-like propagating solution caimot be described with diffusion 
equation. Hence, the cylindrical ion-soliton presented here is not diffusion but a wave 
phenomenon. 

The electron bubble wake is excited by a sub-wavelength impulse of an ultra-short 
driver. In contrast, the ion-wake is excited as the ions undergo sustained interactions 
with the bubble fields within the spatial extent of the wake over several plasma electron 
oscillations. This happens because the electron wake-plasmon oscillations |l2]] have a 
near speed-of-light phase-velocity (jS^ ^ 1) but negligible group-velocity lETI h ^ 

(in the 1-D limit), where, Vth - ^Jk^TJme is the mean electron thermal velocity of the 
background plasma. Therefore a slowly-propagating train of coupled electron plasmons 
is excited in a cold collision-less plasma 11211 . A large difference between phase-velocity 
and the group velocity of the electron oscillations allows sustained field-ion interactions. 
It should be noted that high phase-velocity plasma electron waves are possible only in a 
cold plasma with appropriate density, Uq that allows near speed-of-light propagation of the 
energy sources, jSes - 1 ^ jS(^. Ion-soliton modeled here is assuming a significant difference 
between the phase-velocity and the group velocity of the plasma-electron waves. 

A time symmetric electron wakefield would excite time symmetric ion oscillations 
where the ion velocities average to zero. However, the bubble wake is asymmetric in 
time as the back of the bubble electron compression is a small fraction of the length of 
electron cavitation. The electron oscillations become non-linear at high driver intensities 
as all the interacting electrons are displaced radially, bnJuQ > 1, forming a non-linear 
bubble-shaped electron spatial structure enclosing ions in its cavity. The wakefields 
excited in the bubble are useful for accelerating electrons BUBZlllBl. High intensities also 
lead to fields that can directly drive the plasma electrons to velocities near the speed- 
of-light. This occurs when for a laser pulse Uq > 1 and an electron beam ^ (c^) - ^ 
where Uq is the peak normalized laser vector potential, n^, the peak beam density and 
radius. The radially expelled electrons oscillate radially under the force of the plasma 
ions. These oscillations are excited over plasma electron oscillation timescales, 2na)~} 
{cOpe = where yef^ei^eC is the temporally anharmonic relativistic electron 
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quiver momentum. The normalized quiver momentum of the electrons in the bubble- 
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oscillations is relativistic y±(i± > 1 and the quiver frequency is a)j_ = cOpe |- 

We show fhaf non-linear ultra-relafivisfic elecfron Wakefields inferacfing wifh fhe 
plasma ions lead fo the excitation of a non-linear ion-wake. The non-linear ion-wake 
6n,/no > 1 in Fig ^ and Fig, ^ is excifed over timescales » in the trail of a bubble- 

wake train. By shaping the energy source it can be matched or guided to excite a long train 
of nearly identical plasmons, Figj^ Since it is the electric field Eu,k of a nearly sfationary 
bubble plasmon fhaf excife collective ion-motion we model the ion dynamics in a single 
bubble. Using the single bubble ion dynamics, Figj^we model the ion-wake over the 
whole bubble-train spaiming several hundred plasma skin-depths (c/cope). 

The wake-plasmon energy density {S^k = 0-5{e\Ep\/{meCa)pi;))^ tiq, where Ep is the 
Wakefield amplifude) is continually partitioned befween the field energy and the coherent 
electron quiver kinetic energy. In our model we do not include heavy beam-loading 
of the bubble electron wake. Under heavy beam-loading the bubble field energy is 
efficiently coupled to the kinetic energy of the accelerated beam. In this scenario the 
bubble collapses and the magnitude of fhe ion-wake is smaller. The decoherence of fhe 
ordered elecfron quiver fo random thermal energy, ^nTwk due to the phase-mixing 

| |20| of individual electron trajectories caused by the non-linearities and inhomogeneities 
is further stimulated by the ion motion. The details of the thermalization of the wake 
electrons under ion motion is beyond the scope of fhis paper. It is over these timescales 
upon thermalization that the steepened ion-density expands outwards radially as a non¬ 
linear ion-acoustic wave driven by the electron thermal pressure. The energy transfer 
process observed here is a coupling from fhe non-linear plasma elecfron-mode fo a non¬ 
linear ion-acoustic mode |l23l. We also observe energy coupling to the bow-shock which 
is formed behind the bubble, Figj^ 

The ion-wake is a partially-filled charmel wifh sub-skin-depfh densify-spikes on-axis 
and af fhe bubble-edge locafed af the bubble-radius, Kb IS (6| of several clwpe- The ion 
accumulation in both the density-spikes is many times the background density as shown 
in Figj^c) and Fig|^c). The wake energy is coupled to the ions and to the thermal energy 
IfT^ of the de-cohering wake electrons. The interplay between the electron thermal energy 
and ion kinetic energy then form a cylindrically symmefric non-linear ion-acousfic wave 
which further empties the near-void region of the ion-wake charmel. 
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The time-scale of dissipation of the ion-wake and relaxation of the plasma distribution 
to Vfhlc ~ 0 sets an upper limit on the repetition-rate IITtI of the future plasma colliders. It 
is well known that the ion acoustic wave is damped by collisions and ion-wave Landau 
damping. The ion dynamics also opens important questions on the design the plasma 
container walls in a future high-repetition rate accelerator because if a significant flux of 
ions strikes the wall it can cause significant damage. It was suggested that the wakefield 
energy in the plasma wave could be replenished and sustained Il24ll by a train of energy 
source in order to achieve high repetition rate. However, as shown in this paper due to 
ion motion this is not possible. 

We explore the use of the ion-wake chaimel for positron-beam driven positron accel¬ 
eration in a novel and relevant "suck-in" regime ITSl where the positron beam radius 
Tpi, > cjcope- Such channels are also promising II27111351113^ for exciting electron-wakefields 
with independent transverse and longitudinal field spatial structures. 

The formation of plasma charmels which are vaguely similar to the cylindrical soliton 
shown here but are excited by significantly different processes have been shown previ¬ 
ously. These processes are: using a collimated laser with aimular profile 130111311 . by an 
electrical discharge 1321 . by linear wakes in the self-modulated regime 12911281 and by the 
ponderomotive force of a short-pulse laser-driven linear wake 1261 . In the linear electron 
wake regime the wakefields are symmetric so they average out over the electron oscilla¬ 
tion period and only the second-order ponderomotive force (Fp'‘ ~ -V{Ewk)^) drives the 
ion-motion 1231 . Spatially, in the linear-wake regimes the radial ion excita¬ 

tion is limited to around cjcOpe and the ion-density perturbation is linear in the background 
plasma density, tiq 12611231 . Additionally, in the earlier studies no ion acoustic-waves have 
been examined. 


III. ION-WAKE AS A DRIVEN 
RADIAL ION-ACOUSTIC WAVE 

To develop insight into the ion wake physics, we consider the 1-D simplified dispersion 
relation of the ion-acoustic plane waves. 
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( 1 ) 


where, copi = cOpe yfmjmi and Cs 


1 + {CsjcOpiYk^ 
jTkBTu,k/mi under the collision-less condition, T\ ^ 


and T = 1-1-2//is the adiabatic index with / being the degrees of freedom of the ions. 
During the initial phase of the interaction between the bubble fields and the ions, the ion 
inertia dominates their motion. At times when the fields of the bubble interact with the 
ions, the inertia of ions leads to very small spatial displacement scales and so > oo so 


the term fc(cs/ca„,) » 1 where CsIcOp 


is the Debye wavelength. Thus in this 


phase 0 ) = cOpi and the ion-soliton spike at Rb grows over the plasma frequency timescale, 
2nco~h At initial times the plasma is also cold so Cs - 0. 

When the ions gain significant momentum and start oscillating over larger spatial 
scales in response to the electron dynamics then kAoe ^1- In this phase the acoustic 
wave propagation becomes dispersion-less with co = kCg. However, when the ion density 
perturbation undergoes steepening in a large amplitude acoustic wave forming a spike 
over much smaller spatial scales, k becomes large and dispersion becomes important. 
The ion-acoustic dispersion relation as k increases has higher-order corrections ink, co = 
Csfc - / j ]<Y. Also, as the ions are heated at later time T, increases and modifies the 
sound speed to Cg = ■\lkB{Ti,Te + T,T,•)/?«,. Thus when the ions are heated the sound speed 
lies between and ylTJcYTJmi. 

A. Linearized ion-acoustic wake 


During the time when the oscillations in the electron plasma-wave wake are undergo¬ 
ing phase-mixing the Wakefields still persist and the ion motion is driven by the fields as 
well as the thermal pressure. We show under the linear approximation how an ion-wake 
differs from a freely propagating ion-acoustic wave. It should be noted that as the wake 
electrons phase-mix, the plasma is not isothermal in this case. The linear ion-acoustic 
wave can be obtained by perturbative expansion of ion density, n, and ion fluid velocity 
Vj in the zeroth-order ion fluid continuity equation, HqV • v| -I- = 0. Taking a partial 

pa 

derivative with time, V • + Jp/p = 0- The ion-fluid equation of motion where the 

electron temperature has a spatial gradient is, = eZiE^^k - -T^By-Vr^. 
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FIG. 3; Bubble-wake train behind an ultra-relativistic electron beam under the condition fig(bubble) « fi^bbubble) jibeam with 
beam-plasma parameter as shown in in Fig[^. This figure shows that the bubble which is excited at the earliest time undergoes 
phase-mixing over as several cycles. The intermediate stages of the extent of phase-mixing can be inferred from the bubbles that are 
closer to the beam. Here the beam is located between 170 and 180^ (easily inferred from the radial electric field), (a) electron 
density in 2D cylindrical real-space, (b) corresponding longitudinal electric field profile and (c) corresponding radial-field profile. 
The beam-plasma parameters are the same as in Figlllbut the electron-wake is shown at an earlier time t = ISOaipe- 


Upon substituting the equation of motion in the time-derivative of the linearized conti¬ 


nuity equation, V ■ 


eZ, 


-•wk 


rksT, 


( 1 ) 


tlQ 


Tke " 
m no 


-vrJ+ 


dfi no 


= 0. Thus, a driven ion-acoustic 


wave linearized to the first-order in density perturbation has the form. 


'21 

^ i no 


V- E^;t(r,0 

nii 




wake nti no 


^ I thermal 


( 2 ) 


In this first-order approximation the right-hand side of the eqj^ shows two separate 
timescales of the ion-wake. At earlier times, the first term on the right-hand side domi¬ 
nates. This is the formation phase of the ion-wake where the bubble electron oscillations 
undergo ordered radial motion and the bubble radial electric field excites the inertial re¬ 
sponse of the ions. This leads to the formation of the on-axis and Rb ion density spikes. At 
later times after the phase-mixing between radial oscillators the electrons thermalize and 
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Eit))c(r/ 0 ~ 0. This is the propagation phase where the electron thermal pressure gradient 
drives the cylindrical soliton at Rb radially outwards. 

The PIC simulations show the two phases of the electron dynamics. The coherent 
and phase-mixing phases of the electron oscillations are evident in Figj^a). The electron 
motion in the first oscillation bucket is coherent whereas several oscillations behind the 
driver, the electron begin to de-cohere. This is also evident in the time-evolution of the 
ion soliton in Figj^b) where the radial electric field goes to zero around 200a)~} and the 
ion-soliton is seen moving outwards radially in[^a). 


B. Second-order non-linear ion-acoustic wake 

The linearized ion-acoustic wave equation is however inadequate to describe the prop¬ 
agating solitary density spike at the ion-wake edge, with ion density accumulation many 
times the background density. In the linear regime the homogenous ion-acoustic wave 
equation predicts sinusoidal radial ion oscillations that support the wave. 

When the density in the ion perturbation begins to rise to the order of the background 
density, the electrostatic potential due to charge-separation between the ions and the ther¬ 
mal electrons correspondingly rises. This leads to wave-steepening due to the preferential 
acceleration of ions in the direction of the ion-acoustic wave velocity. When the potential 
of the wave is large enough the ions get trapped and co-propagate with the ion-acoustic 
wave phase velocity, this non-linearity is the basis of the soliton. It should be noted 
that the linearized kinetic theory does not formally incorporate the trapping of particles 
at the wave phase-velocity. In this limit the density perturbation shape is therefore not 
sinusoidal as the co-propagating background ions accumulate and their density pertur¬ 
bation takes the form of an ion-soliton. The co-propagating ion velocity in the soliton can 
therefore exceed the ion-acoustic phase velocity, Vi > Cg and A1 - 1 > 0 where M. = VijCs is 
the Mach number. Therefore, non-linear acoustic waves are in the form of a soliton and 
propagate faster than the ion-acoustic velocity. 

In the second-order, the non-linear ion-density spike n/(r, t) > no propagation is gov¬ 
erned by the Korteweg-de Vries (KdV) equation IlMl which has propagating solutions of 
the form *tl{r- Mcgt) HT^ where *14 is the ion-acoustic waveform, a soliton solution and M 
is the Mach number of the propagating solution. Higher-order contributions to the KdV 


12 


equation have also been considered. However, the standard form of the KdV equation 
assumes an isothermal plasma whereas the bubble-wake phase-mixes into a plasma with 
radial temperature gradient. In a non-isothermal plasma the effect of trapped electrons 
in the ion-soliton have been considered using the Bernstein-Greene-Kruskal model at the 
ion-acoustic velocity IIT^ . 

It is also known that a single ion-soliton under the appropriate conditions can break¬ 
up into multiple solitons leading to a N-soliton solution. This is also a phenomenon we 
observe in the simulations shown in the ion density of the beam-driven case at z = 60-^ 

J COpe 

in Fig|g 

We consider a description of the non-linear cylindrical ion-acoustic waves with a radial 
temperature gradient by retaining the second-order of perfurbafion in ion-densify and 
fluid velocity. We also assume that the background electron trapping does not significantly 
modify the distribution function. We assume that the temperature does not change by 
several orders of magnifude in vicinify of the ion soliton. This assumption is also validated 
by the PIC simulations in Figj^ 



0 2 4 6 8 10 

Radius [c/cOp] 


FIG. 4; Radial profile of the root-mean-square radial electron momentum (proportional to the square-root of the electron 
temperature, ^/T^) at 460 aip} for the beam-driven ion-wake in Figj^ The blue curve shows the root-mean-squared radial electron 
momentum, >')ne{k, r)j r), profile of the wake electrons corresponding to the time in Figj^at 460 Wp}. This 

represents the square-root of the electron temperature, oc The radial gradient of the temperature, ^ T, is thus computed at 
the peak of the soliton (red) and in its vicinity (green). It is interesting to note that = 0. 

\peak 


To obfain fhe KdV equation in cylindrical coordinates with radial temperature gradient 
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time: 


we normalize with respect to the local electron temperature, the radius: ^ 
t = cOpit = t electric field: E = ^E, potential: (p = ion-density perturbation: 

hi = Hi/no, electron-density perturbation: hg = Wg/no, ion-fluid velocity: v = Under 
this normalization the cylindrical coordinate equations transform as: electron Boltzmann 
distribution equation ^ = -hgE - hgCp-^ biTg, ion-fluid continuity equation + 2^ + 
2-^hiV = 0, ion-fluid equation of motion j.v + v^v = E and the Poisson equation = 
4^(fE) = hi — hg. The electric field E is both due to the thermal pressure and the radial 
fields of the wake, £!a,fc -i- But, in the following analysis the propagation of a non¬ 
linear ion-acoustic wave is considered, so we assume that the electron oscillations are 
thermalized and thus the effect of the fields of the wake is negligible, Eu,k —> 0. 

We look for a propagating disturbance of hg, hi, v and E in a stationary background 
plasma with uniform background density no. We consider weakly non-linear ion-acoustic 
wave and expand all the wave quantities in the powers of 6 = A1 - 1. We expand the ion- 
density: hi = 1-I--I-electron density: hg = l+6nf^+6^nf'^ +0{nf^), electric- 
field: E = 6E^^^ -i- + 0{E^^^), electrostatic potential: (p = 6(p^^^ + + 0((p^^^), electron 

temperature: Tg = ^ -i- ^ -i- 0{Tf^) and ion-fluid velocity: v = 6vf^ + + 0{vf^). 

Note that we have assumed that before the electron wake excitation the plasma is cold, 

^ 0. 

We transform to a moving frame of the steepened ion density perturbation using the 
coordinate transform E = - f) and t = Therefore, the radial f-derivative is, 

I = ii = the f-derivative is | = + || = Using this, 

f = + 6 ~^t) and % = = ~l- We renormalize the electric field as, E = 6~^^^E. 

Note that in the moving frame the potential gradient is, E = = so E is a 

more appropriate quantity 

Under the assumption that in the moving-frame the quantities of the disturbance 
change with small 6 = A1 -1, the terms in equations are perturbatively expanded and the 
terms with same powers of 6 are collected. From the 6^ order terms of all the equations 
above, we infer ^ = "tl and = -E^). 

Collecting the 6^ terms from the Boltzmann's equation we obtain E^^^ = --^nf^ + 
From the ion-fluid equation of motion we obtain - j^nf’ = 

jfU + and from the Poisson equation we obtain “ ^P)- Tak¬ 

ing the 6-order terms of the continuity equation and substituting *14 we obtain, *14 + 
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= 0. Neglecting quantities with 6 
times the second-order terms and using the ^'1/ result above, ^ + 2^'W -i- I'Ll-^14 + 
[^nf^ - Using the terms of the Poisson equation we obtain the 

driven Korteweg-de Vries equation in cylindrical coordinates IfTTlI . 

= 'I/ 

1 

It differs from the cartesian-KdV equation by the term ^ arid the driver term 
The cartesian KdV equation can be analytically solved to obtain two classes of solutions: 
(a) self-similar solutions which are shown in IflSB to be Airy functions and (b) soliton 
solutions. A "soliton" is a single isolated pulse which retains its shape as it propagates 
at some velocity, Vsouton- This means that for a solifon-like solution the 14 only depends 
upon the soliton-frame variable, C = ^ - AIcsT and not on space-like ^ and time-like t 
variables separately The solution of fhe carfesian KdV equation in fhis co-moving frame 
is 1i{Q = 3Vs sech2(.^C)llIl. 

The cKdV equation and the driven cKdV equation obtained here caimot be solved 
analytically However, the numerical analysis and experimental verification IIT2l of fhe 
cylindrical-KdV (cKdV) equation show fhaf it supports functions of the form ^ oc sech^(r- 
AiCgt) in the form of a cylindrical ion-solifon. Buf, fhe amplifude of fhe cylindrical solifon 
changes as it propagates. The velocity of the soliton in the cylindrical case is higher than 
in cartesian case HTTI . Since the ion-wake is excited in a non isothermal plasma its velocity 
changes as it is driven. The mean electron temperature reduces because the volume of fhe 
ion-wake expands as fhe soliton moves radially outwards. It is also known to support N- 
soliton solution. We computationally seek the dependence of the non-linear ion-density 
spike on (r - AtCsf)-coordinate. 

It should be noted that such soliton solutions are supported under certain limiting 
condition on the Mach number. At. The strict condition on the existence and stability of 
an ion-solifon arises from fhe requirements on the magnitude of fhe potential required 
for frapping the background ions for a given velocify of the soliton. 

Here we find thaf the speed of ion solifon is approximafely equal fo the ion-acoustic 
speed using the mean temperature. As this is not an isothermal plasma there is no well- 
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defined ion-acoustic speed. So, the ion-acoustic wave is phase-mixed as its phase-velocity 
is spatially varying. However, we confirm from the simulations that this condition is 
strongly satisfied as At > 1. In this scenario the mean electron temperature provides 
an estimate of the sound-speed. The local electron temperature of the ion soliton as 
shown in Figj^ is used to calculate Mach number. At and thus a stability criteria can 
be derived. This problem is represented using the condition on the Sagdeev psuedo- 
potential, 'y((/)) = - (exp((^) - 1 -i- At(At^ - that it has to be a real number. 

This condition is satisfied when At^ - 2(p > 0 therefore, (p < At^/2 and (pmax = A^I2.. 
Using this we find the well-known condition, 

1 < At < 1.6, Vi < 1.6 Cs 
At^ 

O < = 1.28 (4) 

2 

As will be shown later, we find from simulations that the Mach number calculated 
using the mean temperature is well within these bounds and thus the soliton is stable. 

IV. NON-LINEAR ION-WAKE MODEL: 
lON-SOLITON EXCITATION PHASE 

Since the characteristic time of ion-motion is much longer than the electron oscillations, 
the longitudinal field ■ z averages out over the full bubble electron oscillation. So, the 
ions do not gain any net longitudinal momentum. However, atypical radial ion-dynamics 
arise because the radial fields, Eipk ■ f are asymmetric in time as shown in Fig|^and do not 
average to zero, driving an average radial ion-momentum. 

The first stage of the ion-wake formation is controlled by the different time-asymmetric 
phases in Figj^of the bubble radial impulses namely, "suck-in" due to the electron com¬ 
pression in the back of the bubble fduring T^acJc shown in Figj^c), and the "push- 
out" due to the mutual-ion space-charge Coulomb repulsion force during Zcav shown 
in Figj^d). The suck-in force is spatially-periodic at non-linear plasma wavelength, 
Ajvp ~ 2Rb with a duty-cycle D = « 1. In addition to the plasma wake, the propa¬ 

gating energy sources themselves impart impulses such as the laser ponderomotive force 
where = -^V,z|a(r)p and the radial force of the drive beam where 
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Ion dynamics in a bubble 




FIG. 5; Ion dynamics in longitudinally asymmetric phases of the radial forces in an electron bubble, (a) electron density of a 
bubble in 2D cylindrical real-space, (b) longitudinal on-axis profile of the electron density (black), longitudinal field (blue), focusing 
field (red), (c) radial-field profile close to the back of the bubble. This is the focussing "suck-in" phase for the ions, (d) the fields at 
the center of the ion-cavity of the bubble. This is the defocussing "push-out" phase for the ions. 


Fb{^) = -2ne^ni,r. We neglect the driver impulses (below threshold intensity for direct 
non-linear ion excitation [|9 | ||T0 | h because they act on the ions over their sub-wavelength 
short duration. This is unlike the slowly-propagating wake-plasmon bubbles that un¬ 
dergo continual interaction over many plasma periods. The validity of this assumption is 
evident from the laser ion-wake in Figj^ Since the ponderomotive force of a laser driver 
is an outward force for both the electrons and ions, the on-axis density-spike carmot be 
from fhis force. Similarly fhe ion-densify-spike af the radial wake-edge in an electron 
beam driven ion-motion carmot be excited directly by the force of the beam, and is caused 
by the electron wake's radial-edge density compression. 
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The Lagrangian fluid model of the ions in a bubble consists of ion-rings under cylindri¬ 
cal symmetry with mid^ri/df = 'LF^k. The bared-ion region inside the bubble is assumed 
to be a positively charged cylinder under steady-state approximation {Rb > Tbs, back 
of the bubble electron compression radius). The force on the ions from the non-linear 
electron compression bue = ngg » tiq in the back of the bubble and radius tBe, pulls the 
ion rings inward; and within the bubble, the space-charge force of the ions opposes it 
and prevents full collapse. The "suck-in" force on ions is = —Z{lne^nBe^. The 
space-charge force on the ions in the cavity is F^'^ = Zilne^nQTi. The equation of motion is 
mid^Ti/df - - F^‘’‘^hi,ack) = 0 using, coF = Z, ^neHolmu we have. 


d^Yi 




Cu / 

pi I ^Be ^back 'Be 
2 y Hq Tcav 



(5) 


where we have assumed that CTcavIZ^p — 1. Therefore the equilibrium radius where 
the impulses balance out is = rgg D. The ion-rings at r, < collapse inwards 
towards the axis resulting in an on-axis density spike. Whereas the ion-rings at r, > 
move out away from the axis. For m;/Z, > tYip the ion-response is slower but similar. 

When the radially outward moving ion-rings reach beyond Rb, there is excess net 
negative charge of the wake electrons within the bubble-sheath. As a result the radially 
propagating ion rings get trapped and start accumulating just inside the bubble and 
caimot freely move beyond, forming a density compression at Rg. So, the cylindrical ion 
soliton is formed around Rg. This is seen in Figj^and Figj^ where the ion and electron 
density peak at Rg. 

In the non-linear regime Rg » cjcope, the spatial-scale of the ion-wake is over several 
cjcope- This is due to the balance between the opposite radial forces on the electrons at Rg, 
from the driver and the ion cavity 01161. In the laser-driven bubble Fj’J” = ^V,|ao(r)p - 
Fcav = Ine^YiQRB gives Rg ~ {cjWpef Moir)?- when simplified using V,|ao(r)p - fl^/Rg 
and Ye - Uq, Rb - (computationally, ^ 2 cjoope 0)- In the beam-driven bubble 

Ffo(RB) = 2ne^ni,rl/RB - Fcav = 2ne^noRB which gives, Rg ^ where Ap = npnrl 

is the line charge density of the beam. 

We verify the above model using 21D OSIRIS PIC simulations 1251 of the ion-wake in the 
bubble regime by simulating various laser-pulses in cartesian coordinates and electron- 
beams in cylindrical coordinates. We initialize Eulerian specification of the plasma (non- 
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moving window) of pre-ionized Z; = 1 homogeneous density with a density-ramp over 
20 c/ojpe at each of the vacuum-plasma interfaces. We resolve the smallest spatial scale, 
cjoope in the beam case and c/coq in the laser case (laser frequency ooq), with 20 cells. The 
time-step satisfies causality and minimizes numerical dispersion. It is normalized to 
ojp} for beam and to for lasers. We use absorbing boundary conditions tor fields and 
particles. The laser pulse is circularly polarized with radially Gaussian and longitudinally 
polynomial profile Il3^ with flo = 4 (not shown Uq = 1.0 to 40.0), pulse length of 30^, 
matched focal spot-size radius of 40^, and laser frequency coq = lOcOpg. The electron 
beam is initialized with ~ 38,000, rih = 5no (not shown rih = 0.25no to 50no) and spatial 
Gaussian-distribution with cr, = 0.5-^ and m = 1.5-^. 

^ OJpe 


We compare the electron-beam driven ion-wake soliton structure in theory to the 
simulations in Figj^a) and |^a). The observed Rb = 2.45 cjcope (just behind the beam) 
whereas the estimated bubble radius is Rb = = 2.57 cloope {rb - 2.3cr„ the 

assumption rb « Rb is not well satisfied). In Figj^ which is in the propagation-phase, 
the observed ion-wake soliton is located at r - 3.3 c/a)„e at 460 a)~} which is about 1.7—. 

r pc cOp' 

The ion-soliton is excited at an early time around Rb and in the snapshot in Figj^ it has 
propagated outwards. The on-axis density spike drops to a minimum at r^^ ~ OAScIcOpe 
in Figj^whereas the estimated r^ = 0.5 cjcope {fiBelrio ^ 12, £) - 0.1, rgg - 0.5 cjcope). The 
radial ion momentum ^ phase-space in FigjTpb) shows the ions accumulate at the axis 
and the charmel edge at a time corresponding to Figj^b). The ions at the charmel edge are 
seen to have a drift velocity and a thermal spread. The radial electron momentum pr - r 
phase-space in Figj^a) shows that a large density of thermalized electrons are trapped 
within the ion soliton which is confirmed from the density plots in Figj^a). 


In the laser-driven bubble simulations the expected and observed Rb - 4 cjcOpe as 
shown in Figj^a). In Figj^c) the ion-wake soliton is created at r = 4.2 cjcOpe- The 
expected and observed on-axis density-spike radius is r^ = 0.45 c/cope {riBelriQ ^ 8, ^ 0.1, 

rBe — 0.5 c/cOpe). The model for the excitation of this structure of the non-linear wake has 
been verified for a range of laser and beam parameters from quasi-linear to strongly 
non-linear electron wake regime. 
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V. NON-LINEAR ION-WAKE MODEL: 
lON-SOLITON PROPAGATION PHASE 


As described in section IV the electron bubble-wake train fields excite a cylindrical 
ion soliton. Eventually, the electron oscillations phase-mix and thermalize as electron 
thermal energy on the time-scale of about an ion plasma period. The cylindrical soliton 
then propagates radially outwards driven by the temperature gradient as shown in eqj^ 
This is the modified cKdV equation in a non-equilibrium condition when an electron 
temperature gradient also contributes to the force on the ions. 

The charmel-edge density spike, with a form similar fo the cKdV-solution in the r-AiCgt 
frame as shown in Figj^and Figj^to be propagating radially outwards. The propagation 
phase starts around t = 200co~} as the radial electric fields > 0 as shown in Figj^b). 
The propagation phase is evident in Figj^a) where the red curve is the position of the peak 
of the ion-soliton in time. The cylindrical ion-soliton has propagated from rsoiitoni^60a)~}) = 


3.3clcOpe (also seen in Fig,2) to rsoutoni^^OOcOp}) = 4.1c/(Upe which corresponds to an average 
speed of (Vgouton) = 0.0013c. 


We compare the time-averaged soliton speed (Vgouton) to the average speed of sound, 
Cs/c = vl where the average ^ 0.06 from the electron phase-space (not shown). 
This gives Cg - 0.001c (T = 2 for 2D) in agreement with the average soliton velocity. Using 
this time-averaged analysis we see that A1 - 1.3 and so the stability criteria in eqj^is 
satisfied. 

However, as the soliton moves out the volume between the axis and the soliton edge in¬ 
creases thus the electron thermal energy spreads over a larger volume. The re-distribution 
of the thermal energy over a larger volume leads to the reduction in the temperature with 
time. The soliton is not freely propagating but is driven by the radial gradient of the 
electron temperature as shown in eqj^ The soliton speed thus changes in time as shown 
in the red curve of Figj^c). The sound speed also varies with time and it is estimated 
using the temperature at that instant using, Cs{t) = yJkB{Te{t))lmi in terms of the root mean 
squared electron radial momentum. As the plasma is not in thermal equilibrium and its 
temperature varies radially as shown in Figj^ The root-mean-square radial momentum 
represents the temperature at an instant of time and is calculated over radial dimension - 
from the axis to the soliton location: chaimel-{pjp) = p^(r)ng(r)j I'L’^^^neir) oc yu 
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FIG. 6; Time evolution of the cylindrical ion soliton. (a) electron (black) and ion (red) spike radial positions (in terms of cw^}) 
with time and a third-order fit (green) for the position of the ion density-spike of the soliton. (b) radial Wakefields of the electron 
bubble oscillations (in terms of mcCa}pe^~^) at the electron density spike (magenta) and at the ion density spike (blue), (c) radial 
velocity of the ion density spike of the soliton calculated from the third-order fit curve (red). An estimate of the sound speed (green) 
using the mean temperature, between the axis cSc the soliton location (green) and in the vicinity of the soliton peak (blue), in the 
expression Cs = Since the plasma is not isothermal the mean temperature is calculated by averaging the temperature of 

electrons over the indicated spatial region, (d) gradient of the electron temperature at the soliton ion density peak (blue) and in the 
vicinity of the peak (red). The vicinity of the ion density peak of the soliton is defined as shown in Figj^ 


or in the vicinity of the soliton: soliton-{p^^) = p^(r)ne(r)j oc from 

the pr - f phase-space. The instantaneous sound speed, Cs(f) computed with chaimel-(pj^) 
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is in the green curve in FigJ^c) and Cs{t) computed with soliton-{pj^) is in the blue curve 
in Figj^c). The extent of e around the soli ton peak is shown in FigQ 

We compare the curves in (i) red: Vgoutonit) (from the 3rd-order curve-fitting to the 
position of the ion-density peak), (ii) green: Cs{t) from charmel-{pj^(t)) and (hi) blue: Cs{t) 
from soliton-{p^^(f)) in Figj^c). From the comparison it is observed that they are in good 
agreement. It can be seen that the velocity of the soliton estimated using the location of the 
ion-density peak (red) lies between Cs{t) calculated using the average temperature over the 
chaimel (green) which is the upper limit and Cs(f) calculated using the average temperature 
over the soliton (blue) which is the lower limit. Note, that since Vsoiiton{i)ICs = M{t) > 1 the 
blue curve which is lower than the Vsoiuon curve in red, is a better estimate of the sound 
speed in the vicinity of the ion-soliton. The sound-speed with radial momentum averaged 
over the charmel is over-estimated because the thermal electron spread estimation also 
includes the velocity of the electrons drifting with the ion-soliton. 

We also present the radial gradient of the electron temperature, |;Te(r, f) in Fig,^d). It 
is interesting to note from the blue curve in Figj^d) that the temperature gradient at the 
peak of the ion-soliton is zero, = 0. In the vicinity of the soliton peak we see 

that the gradient of the temperature follows the variation in the ion soliton velocity, this 
follows from eqj^ The vicinity of the soliton peak is shown as the green curve overlaid 
on the thermal momentum curve in Figj^ 

In Figj^b) N-soliton formation is observed in the ion-density at around z - 60^. The 
single-ion soliton is seen splitting into several solitons. The N-soliton solution is formed 
because different ion-rings that are driven in the "push-out" phase have a radial position 
dependent defocussing force acting on them, f®‘^(rj) = Zilne^UQri. This is shown in Fig j^d). 
Thus the ion-rings originating at a larger radii from the axis are pushed outwards with 
a force of a higher magnitude and the rings originating at a smaller radii just outside 
are pushed outwards by a smaller force. So, over a longer time the set of ion-rings with 
a higher initial momentum propagate radially outwards at a larger radial velocity. This 
break-up of a single ion-soliton into N-solitons occurs over a longer time-scale because 
the difference in momentum is small compared to the average momentum. 

The thermal momentum, pf at this time is less than one-tenth of the peak wake quiver 
momentum. There are several reasons for the cooling, such as, transfer of the wake energy 
to the ions and the trapped electrons Il34l , escape of the highest energy electrons and un- 
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trapped ions from the channel edge, energy loss to the bow-shock and the re-distribution 
of the energy over an expanding volume. The peak radial ion-momentum is - 0.005 
which shows that not all the radially propagating ions are trapped. The un-trapped free- 
streaming ions at - IcjcOpe can be distinguished from the ions at the charmel-edge in p, - r 
phase-space. 


Radial Phase-spaces 



CO 

^ - 0.05 

CO 

QC 


- 0.10 


0.01 


5 10 15 

Radius [c/tOp] 



-5000 


5000 


4000 

'c 

3000 ^ 

03 



0 


FIG. 7 '. Radial phase-space snapshots of the electron and ion density in Fig|^ (a) electron p, - r radial momentum phase-space 
showing the accumulation of thermalized electrons within the ion-soliton. (b) ion p, - r radial momentum phase-space showing the 

on-axis and ion-wake edge ion accumulations. 
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VI. POSITRON ACCELERATION IN THE 
ION-WAKE CHANNEL 

Positron acceleration using the ion-wake channel is explored in the non-linear suck-in 
regime of positron-beam radii > cjco^e and peak density rip}, > tlq. In this regime 
J'pfc so the on-axis ion density has a radially limited defocussing force. 

It is well known IITSll that in a homogeneous plasma positron beam driven wakes have 
two major problems IIT^ - (i) The plasma electrons collapsing to the axis under the force 
of a positron beam start at different radii. The positron beam radial force driving the 
"suck-in" decreases with the radii. As a result, the electron compression is not optimal, 
(ii) The plasma ions located in the path of the positron beam result in a de-focussing 
force on it. The transport of the positron in a positron-beam driven wake in homogenous 
plasma is not ideal. 

As the ion-wake channel is a practical realization of the hypothesized ideal hollow- 
channel plasma IUBI we examine its excitation by a positron-beam Ill9l . In this section 
we analyze whether the positron-beam driven wake-fields excited in the ion-wake can be 
used for the acceleration and transport of a positron beam. 

We describe the analytical model of the radial electron "suck-in" based excitation of a 
posifron beam wake in the plasma. The equation of motion of the plasma electron rings 
at r from the axis, under the positron beam suck-in force is 

oc ( 6 ) 

where ^ = c(ipbt - z is the space just behind the positron beam with velocity CjSpf, driving 
the collapse. This is a non-linear second-order differential equation of the form, 

r" = /(r,/,e) 

where / is not linear in r. Under the assumption about the positron-beam properties, 
ribpi^) and being constant during the entire interaction of the positron-beam with 

the hollow-chaimel over its full length. So, upon dropping the dependence on ^ the 
equation simplifies fo its special case which has analytical solutions, r" = f{r,r'). The 
solution to this equation is Ifl^ , 
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Vch erf ( ^J\n{rch/r)) = - ^ 


(7) 


where C = ■ Therefore the collapse time-duration is ^coii = -'>'ch We 

note that there is an anomaly that exists in our problem formulation and the solution be¬ 
cause we have not taken into account the space-charge force of the compressing electrons 
as they collapse to the axis and this force balances the suck-in force of the positron beam. 
Under these approximations the collapse time is 


yjn- 


^ ch 


( 8 ) 


(^pe '\j^bp/^0^pb 

Also, note that we have neglected the initial expansion velocity of the chaimel, 
dTch/dt)^^. For optimal compression avoiding phase-mixing, the electron rings should 
collapse over, Tc ^ DAj^pjc where A^p is the non-linear wavelength of the positron-driven 
wake and !D is the duty-cycle of compression phase. So, the optimal chaimel radius is 


- 2 -s/nSD'^^rpi,. The scaling of the with positron beam parameters is show in 


ch 


Using 2-|D pic simulations in a moving window we study the positron beam driven 
Wakefields in cylindrical geometry. We compare positron acceleration in an ideal (Heav¬ 
iside density function, noH(r - rch)) and an ion-wake charmels (on-axis and channel-edge 
density-spike, channel minimum density of O.lno) with Vch = 2.5 cjcOpi,- For non-linear 


wake parameters Vpi, = l.AcjcOpe, nph = 1.3no and - l.AcjcOpe = 0.25). Fig 

shows that the peak on-axis accelerating field is 0.4 m.eCC 0 pee~^ for an ideal channel and 
0.2 meCa)pee~^ for the ion-wake channel. Figalso shows that the focussing potential (nor¬ 
malized to 27.6 m.f,c^e~^) is similar and overall focussing in both cases. However, in the 
ion-channel the radial field is defocussing around the on-axis ion-spike. 

Thus the ion-wake channel is useful for accelerating and transporting positrons despite 
the lower accelerating fields in comparison with the ideal channel. We also note that the 
on-axis density spike has detrimental effect on the focussing fields near the axis. Ideal 
channels of a few c/oipg are technologically challenging whereas the ion-wake channel of 
radius Ych ^ Pb is formed behind every bubble-wake. 
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Positron-beam driven wakefields in an ion-wake channel 



FIG. 8; Radial profile of wakefields excited by a positron-beam in ideal-channel versus an ion-wake channel. The radial 
profile of the normalized electron density (black) in an ion-wake channel (normalized to the maximum electron compression) at 
longitudinal location of the peak accelerating wakefield {rpi, = IScjwpe, yph = 38000, = 1.3no)- The radial profile of the 

accelerating-wakefield and normalized focussing-wakefield potential (radial field integrated from the edge of the box to a radius). 


VII. CONCLUSION 


In conclusion, using theory and PIC simulations we have shown the dynamics of 
the formation and evolution of a non-linear ion-wake excited by the well-characterized 
electron bubble-wakefields Il3l-|l8]]. We have shown that the non-linear ion-wake has the 
characteristic cylindrical ion-soliton solution and evolves to an N-soliton solution over 
longer time as described by the cKdV equation. Thus over the period of persistence of the 
ion-soliton a second electron bunch carmot be accelerated in the plasma. This establishes 
an upper limit on the repetition rate of the plasma collider. We have also shown the 
feasibility of using the ion-wake chaimel for positron acceleration in an experimentally 
relevant parameter regime. 
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